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Abstract. Given interpolation points P\,P2,... ,Pn in the plane, it is known that there 
does not exist an interpolating curve with minimal bending energy, unless the given points 
' lie sequentially along a line. We say than an interpolating curve is admissable if each piece, 

connecting two consecutive points Pi and Pi+i, is an s-curve, where an s-curve is a planar 
curve which first turns at most 180° in one direction and then turns at most 180° in the 
1-^^ , opposite direction. Our main result is that among all admissable interpolating curves there 

d ' exists a curve with minimal bending energy. We also prove, in a very constructive manner, the 

existence of an s-curve, with minimal bending energy, which connects two given unit tangent 
vectors. 
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^ ■ 1. Introduction 

O _ 

^ ' Given functional data {xi^Ui) G R^, it is well known that the natural cubic spline in- 

terpolant y = f{x) (satisfying i/i = f{xi)) minimizes j^{f" {x)Y dx among all smooth 

. functions which interpolate the given data. The functional J^(/"(x))^ (ix is often viewed 

as a simple approximation of the curve's bending energy jj" ^(■s)'^ ds, where s denotes ar- 
clength and k denotes signed curvature, and it is natural to ask what would happen if one 
tried to minimize the bending energy amongst all smooth interpolating curves. Unfortu- 
nately, such optimal curves do not exist except in the trivial case when the interpolation 
points lie sequentially along a line. Apparantly, this was first observed by Birkhoff and 
de Boor [1], along with Birkhoff, Burchard and Thomas [2]. This lack of existence can be 
understood as a consequence of the effect that scaling has on bending energy: the bending 
energy of a curve scaled by a factor A, equals j times the original bending energy. As a 
result, it is always possible to construct smooth interpolating curves with arbitrarily small 
bending energy. For example, one can do this with a series of 'concentric circles' having 
slowly varying, but very large, radii. 



J3 



1991 Mathematics Subject Classification. 41A15; 65D17, 41A05. 

Key words and phrases, spline, nonlinear spline, elastica, bending energy. 



Typeset by AmS-TW- 



2 



ALBERT BORBELY & MICHAEL J. JOHNSON 



Subsequent attention was directed towards interpolating curves whose bending energy 
is locally minimal (ie minimal among all nearby interpolating curves). It was reported in 
[2], and mentioned in [1], that if an interpolating curve F has a locally minimal bending 
energy, then each segment of F, connecting two consecutive interpolation points, will be a 
segment of 'rectangular elastica', meaning a planar curve whose signed curvature k satisfies 
the differential equation 2^ + — Q. (Rectangular elastica was first described by James 
Bernoulli (1694) and resides as one specie among the nine species of elastica identified 
by Euler (1750), sec [4].) Using a variational calculus and physical reasoning, Lee and 
Forsyth [6] (see also [3]) have confirmed that each segment of F is indeed a segment of 
rectangular elastica, and have moreover shown that the signed curvature of F is continuous 
throughout the curve and vanishes at the endpoints. Unfortunately, for a given sequence 
of interpolation points, interpolating curves with locally minimal bending energy do not 
always exist and this constitutes a significant defficiency in the theoretical foundation of 
this interpolation method. 

Rather than seeking an interpolating curve with a locally minimimal bending energy, an 
alternate approach is to define a restricted class of 'admissable' interpolating curves and 
then seek a curve with minimal bending energy in the restricted class. Birkhoff proposed a 
restriction on length and conjectured that among all smooth interpolating curves of length 
at most Lq, Lq being a prescribed constant, there exists a curve with minimal bending 
energy. This conjecture was eventually proved by Jerome [5]. 

Rather than a restriction on length, we propose a restriction on shape. The motivation 
for our restriction comes from the fact that if a smooth interpolating curve F has a locally 
minimal bending energy, then it can be shown that each segment of F, connecting two 
consecutive interpolation points, is what we call an s-curve. In brief, an s-curve is a curve 
which first turns in one direction (either counter-clockwise or clockwise) at most 180° and 
then turns in the opposite direction at most 180°. Given a sequence of points Pi, P2, ■ ■ ■ , -Pn 
in M?, an interpolating curve F is deemed admissable if each piece of F, connecting two 
consecutive interpolation points Pi and Pi+i, is an s-curve. The family of all admissable 
interpolating curves is denoted A{Pi, P2, . . . , Pn), and our main result is the following. 

Theorem 1.1. Given any sequence of points Pi, P2, . . . , Pn in M^, the family of admissable 
interpolating curves A{Pi,P2, ■ ■ ■ ,Pn) contains a curve with minimal bending energy. 

An essential sub-problem which arises in the proof of Theorem 1.1 is that of proving 

the existence of an s-curve, with minimal bending energy, which connects two given unit 
tangent vectors. In addition to facilitating our proof of Theorem 1.1, we anticipate that this 
sub-problem sits at the core of any numerical algorithm for solving the general problem, 
and with this in mind, we present a thorough analysis of the sub-problem along with a 
constructive solution. 

An outline of the sequel is a follows. In section 2, we explain the notation used through- 
out and develop some basic formulae and properties of rectangular elastica. A curve which 
turns at most 180° in one direction is called a c-curve and in section 3, we show the exis- 
tence of an optimal c-curve connecting a unit tangent vector to a line as well as connecting 
two unit tangent vectors. In section 4, the uniqueness, or lack thereof, of the optimal c- 
curves found in section 3 is treated. The important sub-problem mentioned above, namely 
the existence of an optimal s-curve connecting two unit tangent vectors, is primarily solved 
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in section 5, except that one particular case (where the optimal s-curve turns out to be a 
unique c-curve) is treated in section 6. Finally, in section 7, we prove Theorem 1.1. 



2. Notation 

A curve is a difFerentiable function / : [a,6] — > C whose derivative /' is absolutely- 
continuous and non-zero. The length of / is len(/) = £ \f'{t)\ dt. With L = len(/), let 

the variables t G [a, b] and s e [0, L] be related by s = \f'{r) \ dr and define F : [0, L] — >■ C 
by F{s) = f{t). It can be shown that F is a curve (ie F' is absolutely continuous) satisfying 
= 1. The curve F is called the unit speed curve described by / and is denoted 
[/]. Two curves / and g are said to be equivalent, written / = if [/] = [g]. Since F' 
is absolutely continuous and = 1, it follows that there exists an absolutely continuous 
function r : [0, L] — t- R, unique modulo an additive constant in 27rZ, such that F' = e*"^. 
We refer to r at the direction angle of F, while the derivative of r, denoted is called 
the signed curvature of F. Since r is absolutely continuous, it follows that k is Lebesgue 
integrable (see [7, pp. 108-112]). The turning angle of /, denoted A(/), is defined by 

A(/) = A(F) := k{s) ds. Note that the magnitude of the turning angle is bounded by 
the Li-norm oi k. If k > (resp. k < 0) almost everywhere in [0, L], then |A(F)| = 
and / is called a left-curve (resp. right-curve) . A c-curve is a left-curve or a right-curve 
whose turning angle has magnitude at most n. A u-turn is a c-curve whose turning angle 
has magnitude tt. 

Given the signed curvature k of F and its initial position and direction, we can recover 

F as follows: 

Step 1. Define r(s) = tq + /J* K{r) dr, s e [0, L], where tq = arg(F'(0)). 
Step 2. F(s) = F(0) + e^^^'^) dr. 

This reconstruction can be used to show that two curves are 'close'. To see this, suppose 
Fi is a unit speed curve having the same length and initial position and direction as 
F. It follows from step 1, that |T(r) — Ti(r)| < \\k, — kiH^^ and then using the Lipschitz 
continuity of the function r i-> e*^ in step 2, we obtain 







(2.1) \F{s)-Fi{s)\< I e^^^")-e*^i(^) dr< / \T{r) - Ti{r)\ dt < s\\k - m 







iLi 



Whereas the Li-norm of n is necessarily finite, the L2-norm may or may not be finite. 
When it is finite, we say that / has finite bending energy, where the bending energy of 
/, denoted ||/|| , is essentially the square of the L2-norm of k: 



\Ff:=lJ^\K{s)\'ds. 



The constant j has been inserted for later convenience. 

A unit tangent vector is an ordered pair of complex numbers u = {ui,U2) G such 
that 1^2 1 = 1 and can be visualized as the directed line segment, of unit length, having 
position (or base-point) pos{u) — ui and direction dir{u) — U2- For any t e [a, 6], the unit 
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tangent vector to / at t is f{t) — (/(t), |) The unit tangent vectors u = f{a) 
and V = f{b) are called, respectively, the initial and terminal unit tangent vectors of /, 
and we say that / connects u to v. We also say that / connects w to £ if £ is the line 
through f{b) which is parallel to f'{b). If (7 is a curve whose initial unit tangent vector 
equals the terminal unit tangent vector of /, then [/] can be extended by [g] obtaining 
a unit speed curve, denoted f \-\ g, whose initial and terminal unit tangent vectors equal 
those of / and g, respectively, and whose bending energy satisfies ||/U^||^ = + 

A similarity transformation is a mapping T : C — > C of the form T{z) = ciz + C2 
or T{z) = c\z -\- C2, where Ci,C2 are complex constants, Ci ^ 0. The first form preserves 
the orientation (left or right) of a curve while the second form reverses it. The dilation 
factor is A = |ci|, and the effect on a curve / is, as expected, len(T of) — Alen(/), 
||To/||^ = i||/f and |A(To/)| = |A(/)|. If a curve g is equivalent to T o /, then we 
say that g is similar to /; in case A = 1, T is called a congruency transformation 
and we say that g is congruent to /. Furthermore, we say that g is directly similar (or 
congruent) to / if T is orientation preserving. 

The curves constructed in this article are formed by line segments (denoted [A, B]) and 
various segments of rectangular elastica. For the latter, we employ the parameterization 

sin^ t 

E{t) — sint + ii{t)^ where ^{t) is defined by — = — =, ^(0) = (see Figure 6.1a). 

Vl + sin^t 

Since ^ is even and 7r-periodic, it follows that ^ is odd and satisfies + ^) = d-\- ^{f?), 
where d := ^(tt). Since the sine function is odd and 27r-periodic, we conclude that E{t) 
is odd and satisfies E{t + 2tv) = i2d + E{t). We use the notation £"[0,6] to denote the 
sub-curve E{t), t e [a, 6], and any curve which is similar to E]^a,b] is called a segment of 
rectangular elastica. For later reference, we mention the following. 

\E'{t) \ = = cost^/l + sin^t + ^sin^t, K(t) = 2sint, 

Vl + sin^t 1^(^)1 



cos h 



A{E[o,b])^ n{t)\E'{t)\dt = 2 cos-^ 
\\Eia,b] f = \f^ <tf\E'{t) I dt = ab) - ^(a) 



TT 

2' 



For to e (0,7r], the segment -E[o,to] P^^ys an important role in the sequel. In the following 
lemma, we establish a connection between the turning angle of E[q and the value of 

Lemma 2.1. Let to & {0,7t] and put tq = A{E^q^i^]). Then ^{to) = - / Vsinrdr. 

2 Jq 

Proof. Fix to e (0,7r] and put r = A{E[o,t]) = 2cos-i ^ - f , t e [0,to]- Then 
^ = K{t)\E'{t)\, and since e'^ = E'{t)/\E'{t)\, we have sinr = ^E'{t)/\E'{t)\ = sin^t, 
which implies V sin r = sint. Hence, 

r^o sin^ t 1 1 

e(to) = / , d,t=- sm{t)K{t)\E'{t)\ dt=- V^dr. 

Jo VI + sin"' t Jo '''Jo 
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3. Existence of optimal c- curves 



Given a unit tangent vector u and a line £, let Ci{u,£) denote the set of left c-curves 
which connect u to i. In this section, we consider the problem of finding a curve in Ci (w, £) 
which has minimal bending energy. We first consider Ci{uo,£d), where uq = E{0) and 
•^d = {-2 e C : $52; = d}. We will show that -E[o,7r] has minimal bending energy in Ci{uQ,£d)- 

Note that ||£^[o,7r]||^ — ^{'^) — d and by Lemma 2.1, we have d — ^{tt) = \ Jq \/ sin r dr. 

Let / G Ci{uo,£d), put L = lcn(/), and let F = [f] denote the unit speed curve 
described by /. Let r and k be the direction angle and signed curvature of F, respectively, 
and note that k{s) ds = n since the turning angle in Ci{u,£d) is tt. Futhermore, since 
F originates at and terminates on £d, we have d = '^F{L) = sinr ds. 

1 1 2 

Lemma 3.1. If k, is continuous and positive, then \\F\\ > d. 

Proof. We adopt the viewpoint that r e [0, tt] and s e [0, L] are variables related by 
T = Jq K{r) dr. The assumptions on k ensTire that r and s are increasing func- 
tions of one another. Noting that ^ = we observe that ||-F||^ = \ f^^ nis^ ds = 
1 lo '^(^)^ ds = lJ^K{s)dT. Similarly, since ^ = we have Jq ^ dr = 



Jo 



siiiT^ dr = Jq sinrds = d. Now, 



d^ = 



- / V sin T dr 
2 Jo 



-\/sinT \/ k{s) 



dr 



< 



o 2 2 

by the Cauchy-Schwarz inequality. Hence c? > , and therefore ||F|| > d. □ 

1 1 2 

Returning now to the general case, suppose, by way of contradiction, that ||F|| < d. 
Then k is a nonnegative square integrable function satisfying ds — tt. It follows 

that for every e > 0, there exists a positive continuous function Kg : [0, L] [0, 00) such 



that Ke{s) ds = TT and 



< e. Let Fe be the unit speed curve having signed 



curvature Kg and initial unit tangent vector uq. Since F^ has turning angle tt, it follows 
that the terminal unit tangent vector of F^ is parallel to £d, but there is no guarantee 
that the terminal point = Fg{L) lies on £d. We repair this by multiplying Fg with the 

2 111 1 1 2 

positive scalar Cg = df^Zg obtaining the curve CgFg G Ci{uQ,£d) with HcgFgH = — . 



Since \\k 



0, it follows that llFgH as £ ^ 0. And since the Li-norm of 



K — Kgis bounded by a constant multiple of its L2-norm, it follows from (2.1) that Cg — > 1 
as £ — > 0. Hence ||ceFe|| < d when e > is sufficiently small, contradicting Lemma 3.1. 
Therefore, ||F|| > d. We have thus proved the following. 

Theorem 3.2. Let uq = E{0) and £d = {z e C : = d}. 

bending energy in Ci{uo,£d). 



Then -E'[o,7r] has minimal 
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Corollary 3.3. Let ti e (0, tt) and put ui — E{ti). Then -Ej^^^tt] has minimal bending 
energy in Ci{ui,£d). 

Proof. If f E Ci{ui,id) has bending energy less than E^^^^^^, then -Ep,*!] U / will be a curve 
in Ci{uo,£(i) with bending energy less than d, contradicting Theorem 3.2. □ 

Corollary 3.4. Let < ti < t2 < tt and put u\ = E{t\), U2 = E{t2). Then -Ej^^^^^] has 
minimal bending energy in Ci{ui,U2)- 

Proof. If / e Ci{ui,U2) has bending energy less than E[ti^t2], then E[o,t,] U / U E[t^^^j will 
be a curve in Ci{uo,£d) with bending energy less than d, contradicting Theorem 3.2. □ 

Definition 3.5. Let £ be a line and u a unit tangent vector, with pos{u) ^ £, and assume 
Ci{u,£) is nonempty. Let S e (0, tt] be the common turning angle in Ci{u,£) and let 
ti e [0,7r) be such that A{E[i^^.j^]) = 6. There exists a unique similarity transformation T 
such that T o E[i^^^-^ belongs to Ci{u,£). We define l{u,£) = T o E^f^^^y In other words, 
l{u,£) is the unique curve in Ci{u,£) which is similar to E^^^^^^y 

Theorem 3.6. Let £ be a line and u a unit tangent vector, with pos{u) ^ £, and assume 
Ci{u,l) is nonempty. Then l{u,£) has minimal bending energy in Ci{u,£). Moreover, 
if S E (0, tt] denotes the common turning angle in Ci{u,£) and p denotes the orthogonal 
distance from pos{u) to £, then 



p 



V sin T dr 



Proof. Let ti and T be as in Definition 3.5 and ui and £d as in Corollary 3.3. Note that T 
maps Ci (ui , £d) onto Ci (u, £) and a curve / has minimal bending energy in Ci{ui, £d) if and 
only if T{f) has minimal bending energy in Ci{u, £). It therefore follows from Corollary 3.3 
that l{u,£) — To -Eitj^Tr] has minimal bending energy in Ci{u,£). In order to compute the 

2 1 2 

bending energy of l{u, £), recall that £) \\ = i^||i?[t^^7r] || , where A is the dilation factor 
in T. Since the orthogonal distance from pos{ui) to £d is d — ^(ti), it follows that A = 
p/{d—^{ti)). Andsince = (i— ^(ti), we have ||/('U, ||^ = ^(d— ^(ti))^. By Lemma 

2.1, ^{ti) — Jq ^ V sin T dr and hence d — ^{ti) — ^ J^^g V sin r dr — ^ Jq V sin r dr. □ 

Remark 3. 7. The definitions and results for right c-curves are analogous to those of left 
c-curves. In brief, we denote the set of right c-curves connecting w to £ by Cr{u,£), and 
r(«,£) is defined the same as l{u,£) except that 6 denotes the magnitude of the common 
turning angle in Cr{u,£) (right curves have a negative turning angle). Theorem 3.6 then 
holds with Cr{u,£) and r{u,£) in place of Ci{u,£) and l{u,£), respectively. 
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4. Uniqueness of optimal c-curves 

Having settled the question of existence of an optimal curve in Ci{u,i) we now address 
uniqueness. As with existence we start with C/(tto,-^d), where uq and £d are as in Theorem 
3.2. 

Theorem 4.1. For i = 1,2, let Fi : [0,Lj] — >■ C 6e a unit speed curve in Ci{uo,£d) such 
that \\Fi\\ = d and assume that Fi does not begin or end with a line segment. Then 
Fi = F2. 

Our proof of this employs the following technical result, which is left as a simple exercise 
in differential calculus. 

Lemma 4.2. Let Ai, A2 > and define H : (0, 1) (0, 00) by H{u) = — -\ —. Then 

H 1 — A* 

H has a unique minimum at = Ai/(Ai + A2), where H^hq) = (Ai + A2)^. 

Proof of Theorem 4-i- Let Ti and Ki be the direction angle and signed curvature of Fi, 
respectively. Since Fi does not begin or end with a line segment, we have < Ti{s) < tt 
for all s G (0, Li), and it follows that Fi can be reparameterized as t t— > 9i{t) +it, t & [0, d], 
where gi is continous on [0,d] and continuously differentiable on (0, d). Fix 7 G (0, tt) 
and let t G (0, tt) be such that A{E[Q^t]) = 7- Let Sj G (0, L^) be such that ri{si) = 7, 
and put Vi — Fi{si) and ti = '^Fi{si). We claim that ti = ^{t) = ^2- Noting that the 
turning angle in Ci{vi, id) is tt — 7 and the orthogonal distance from pos{vi) to £d is d — ti, 
and since -F^s.^l] belongs to Ci{vi,£d), we obtain from Theorem 3.6 and Lemma 2.1 that 
ll-^i[si, Li] 11^ ^ ~ ^(^))^- By a similar argument (using right c-curves) we obtain 

\\Fi[o,s,]f>j:m''- Therefore, 

d=\\Fif = ||F,[o,.,]||'+ > + ^((i-e(t))'. 

With Ai = ^{t), A2 = d — ^{t), 11 = ti/d, and with H{ii) as in Lemma 4.2, we can express 
the above inequality as (i > ^H[ii), or equivalently, > H{fx). By Lemma 4.2, H has 
a unique minimum at fio = ^{t)/d where H{po) = d"^. But since d'^ > H{fi), it must be 
the case that n — Hq; therefore ti = $,(t) as claimed. In terms of the functions gi and g2, 
we have proved that if gi{ti) — cotj — 5'2(^2), then ti — ^{t) — t2- Since, for i= 1,2, g[ 
is continuous and decreasing on (0,(i), with range (—00,00), we conclude that g'^ = g'2 on 
(0, d). Since ^i(O) = = 5f2(0), we have gi = g2 on [0, d]. From this we conclude that Fi 
and F2 are equivalent, but since both are unit speed curves, they must be equal. □ 

As an immediate corollary, we have the following. 

Corollary 4.3. /// G Ci{uo,id) has minimal bending energy, then f contains a subcurve 
which is equivalent to c + E^q^^^] for some constant c > 0. 

Imitating the proof of Corollary 3.3 and Theorem 3.6, one easily obtains the following. 
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Corollary. Let i be a line and u a unit tangent vector, with pos{u) I, and assume 
Ci{u,i) is nonempty. Let d G (0, tt] denote the common turning angle in Ci{u,£) and let 
f G Ci{u,£) have minimal bending energy. 

(i) If d = TT, then f contains a subcurve which is congruent to l{u,l). 

(a) If 5 < TV, then either f = l{u,£) or f = l{u,l) U [A, S] for some line segment [A,B]. 

We have seen in Corollary 3.4 that E^ti,t2\ has minimal bending energy in Ci{ui,U2). 
Using the same technique as above, one easily obtains the following. 

Theorem 4.4. Let < ti < ^2 < tt, put ui = E{ti), U2 = E{t2) and assume f G 
Ci{ui,U2) has minimal bending energy in Ci{ui,U2). If f is not equivalent to E^^^^-i-^-^, 
then [^1,^2] = [0, tt] and a sub-curve of f is congruent to -E'p.Tr] (i^ f equivalent to 
[0, c] U (c + -E'[o,7r]) U [c + zd, id] for some c> Q). 



An s-curve is either a c-curve (considered a degenerate s-curve) or a curve of the form 
/ = /i U /2, where /i and /2 are c-curves which turn in opposite directions. Let u and v 
be two unit tangent vectors and let S{u, v) denote the set of all s-curves which connect u 
to V. In this section and the next, we will prove the following. 

Theorem 5.1. Let u and v be two unit tangent vectors with distinct positions. If S(u,v) 
is nonempty, then there exists a curve in S{u,v) with minimal bending energy. 

In addition to proving existence, our proof of Theorem 5.1 will actually describe all 
optimal curves in S{u,v). Expecting that the numerical problem of finding an optimal 
curve in S{u, v) lies at the heart of future algorithms, we have structured our proof so that 
it easily translates into a numerical algorithm. 

To begin, let u and v be two unit tangent vectors with distinct positions. By applying 
a similarity transformation, if necessary, and possibly a direction reversal (ie S{—v^ —u) in 
place of S{u, v)), we can assume without loss of generality that u = (0, e*") and v = (1, e*^), 
where a G [0, tt] and |/?| < a (see Figure 5.1). 



We leave it to the reader, as a worthwhile exercise, to verify that S{u, v) is non-empty if 
and only if a < tt and (3 > a — Tc. With that in mind, we proceed assuming that a G [0, tt), 
< a and /3 > ct — tt. 



5. Optimal s-curves, part I 
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If q; = 0, then ^ = as well and the line segment [0, 1] is the unique curve (modulo 
equivalence) in S{u,v) having minimal bending energy. Having dispensed with the trivial 
case, we assume henceforth that a > 0. 

Our proof of existence will show that there exists an optimal curve in S{u, v) having 
one of the following two forms. 

Definition 5.2. A curve / is of 

(i) first form if there exist — tt < to < t < tt such that / is directly similar to -E'[to,t]7 

(ii) second form if there exists c > and t E [0, tt] such that / is directly similar to 

£;[_,,o]U[o,c]u(c + £;[o,t]). 

Note that curves of first form do not contain u-turns, while curves of second form do. 
While studying right- left s-curves in S{u,v), the following quantities will gradually take 
on significance, but for easy reference we gather and define them here. In relation to a 
generic right- left s-curve / in S{u, v), the angle 7 is illustrated in Fig. 5.3, and the set of 
all possible angles 7 is denoted F. 

r [a - TT, /3] if /3 < 
Definition 5.3. For 7 in F := < , we define the following: 

{ [a — n,0) if ^ > 

yi ■= Viil) ■= 2 J Vsin rdr (bending energy of -E[o,ti]) 
y2 ■■= 2/2(7) 2 J Vsin rdr (bending energy of -Ep^t^]) 



1 o 2 

G{j) := : {yi + 1/2) (lower bound on ||/|| ) 

— sm7 

sin 7 



(7(7) := cos 7 H {\/sm{a — 7) + \/sm( f3 — 7)) (signed distance) 

yi + 2/2 

A(7) := — ^^^^ (dilation factor) 

yi + 2/2 



Note that, by Lemma 2.1 (see Fig. 2.1), yi and j/2 can also be expressed as yi = ^(ti) = 
llEpQ^^^jll and y2 = ^(t2) = ll-^fo.ta] H' "^^^^^ *i'^2 e [0,7r] are determined by A{E[o^t^]) = 
a - 7 and A{E[o^t^]) = ^ - 7. 

We mention further that G{'y) (sec Theorem 5.6) is a lower bound on the bending 
energy of our generic curve /, and (7(7) is a signed distance, which is illustrated in Fig. 
5.4. Regarding A(7), we mention that the curves r{u,i) and l{£, v), shown in Fig. 5.4, are 
similar to £"[0,*!] and -Eicta]? respectively, with common dilation factor A(7). The crucial 
identity relating (^(7), c'"(7) and A(7) is given in Lemma 5.10. 

Our constructive proof that S{u, v) contains an optimal curve is broken into several 
cases which depend on a and /3. To help the reader track these cases, we give here a short 
description of each case and where in this section or the next it is treated. 
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Summary 5.4. We assume a e (0,7r), < a and ^ > o; — tt. 

(a) ^ = a — TT. This case is treated just above Remark 5.5 and results in an optimal curve 
of second form. 

(b) /5>0or(Q; — 7r</3<0 and cr(/3) > 0). It is shown in Lemma 5.11 that the function 
G has a minimum value Grain-, and in Corollary 5.12 {vi) it is shown that Grain equals the 
minimum bending energy in S{u^ v). Each 7 G F, where G is minimized, gives rise to an 
optimal curve in S{u,v), but the form of the optimal curve depends on whether or not 7 
is the left endpoint of F. If G is minimized at the left endpoint 7 = a — tt, then it is shown 
in Corollary 5.12 {in), that the curve fa-ir, which is of second form, is an optimal curve 
in S{u,v). If G is minimized at any other point 7 > a — tt, then it is shown in Corollary 
5.12 (iv), (f ) that the curve f^, which is of first form, is an optimal curve in S{u, v). 

(c) a — TT < /3 < and cr(^) < 0. In Theorem 6.2, it is shown that the unique curve 
(modulo equivalence) in S{u, v) having minimal bending energy is a c-curve of first form. 

The case a — tt < /3 < and a(/5) = 0, which is included in both cases (b) and (c) 
above, serves as a bridge between these two cases; the first half of section 6, culminating in 
Corollary 6.8, is dedicated to showing that the function G is uniquely minimized at 7 = ^ 
when a — TT < /3 < and a{/3) = 0. 

Our analysis employs an initial partitioning S{u, v) = S'i^{u, v)USri{u, v), where S[^{u, v) 
(which is non-empty if and only if a — tt < /? < 0) denotes the set of all non-degenerate 
left-right s-curves in S{u, v) and Sri{u, v) denotes the set of all right-left s-curves in S{u, v). 
The set Sri{u, v) will be further partitioned into subsets Sj{u, v), 7 G F, and for each 7 e F 
a distinguished subset s*{u, v) C s^{u, v) will arise. The case 7 = < is exceptional, so 
we address it first. 

When P < 0, Cr{u,v) is non-empty and we define s*p{u,v) := Sfj{u,v) := Cr{u,v). Let 
ij3 denote the line through the point 1 with direction e'^^. Since Cr{u,v) C Cr{u,i/3), it 

2 2 

follows that 11/11 > ||r(t(,, II for all / G Cr{u,v). Since the turning angle in r{u,£f3) is 
— {a — /3) and the orthogonal distance from to £(3 is — sin /3, it follows from Theorem 3.6 
and Remark 3.7 that ||r('U,£/3)|| = G{(3) and that r{u,£i3) is directly similar to E[_f^^Qy 
It can be shown, in a manner similar to the proof of Theorem 3.6, that r{u, Ip) is directly 
congruent to A(/3)-E'[_t^ o]- Furthermore, if Pi denotes the terminal point of r(tt, Ip), then it 
can be shown that the signed distance, in the direction e*'^, from Pi to 1 is cr(/3). We omit 
the details of this calculation since a very similar calculation will be performed below. The 
case 13 = a — TT is particularly easy to handle since F = {a — tt}, rather than an interval. 
Assuming 13 = a — tt (which implies a > 7r/2), we note that the turning angle in r{u,£p) 
equals —n and therefore S{u,v) = Cr{u,v). One easily verifies that (t{P) = — coso; > 
and it follows from the results of section 4 that r{u, £13) U [Pi, 1] is the unique curve, modulo 
equivalence and elongation of u-turns, in S{u, v) with minimal bending energy. Moreover, 
r('U, £/3) U [Pi, 1] is of second form with c = a{l3)/X{j3) and t = 0. 

Remark 5.5. If a curve / G S{u,v) contains a u-turn (eg. the curve r{u,£j3) U [Pi, 1] 
above), then it is always possible to elongate the u-turn by inserting a pair of congruent 
line segments before and after the u-turn. Although longer, the resulting curve still belongs 
to S{u,v) and has the same bending energy as /. 

Having treated both the trivial case a = and the extreme case (3 — a — n, we proceed 
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assuming that 



(5.1) 



a G (0, tt), |/?| < ck, (3 > a — it. 



Note that F is now an interval with right end point (3o := min{/3, 0}. 
Proposition 5.6. // /3 < 0, then ||/f > G'(/3) for all f G S'i^{u,v). 

Proof. Let / : [0, L] — )• C be a non-degenerate left-right s-curve in S{u, v). There exist < 
So < S2 < L such that A = /(sq) is an inflection point (ie /[o,so] is a left c-curve and /[sp,L] 
is a right c-curve) and C = /(S2) lies on the line through the point with direction e*" 
(see Figure 5.2). It can be shown that arg/(so) < arg/'(so) while arg/(s2) > arg/'(s2), 
and it follows by continuity that there exists si G (sq, S2) such that arg/(si) = arg/'(si). 
This equality implies that g := [0, S] U f[s^,L] is a right c-curve, where B = f{si). With 
cti = arg/(si) and ui = (0,6*"^), we note that g G Criui^ip) with bending energy 

2 2 2 

WqW = ll/[s(i,L]|| < 11/11 • Since si < S2, it follows that ai > a, and therefore 



^>\\9f>\\r{ui,ip)f 



1 



1 



sin/? \2 Jq 



ai-/3 



V sin T dt 



> 



1 



1 



— sin/3 \ 2 Jq 



a-fi 



V sin r dt 



\r{ujp)f = G{(5). 



□ 

Let / be a non-degenerate right-left s-curve in S{u,v). Then / has a well-deflned 
inflection line I and inflection direction 7, and it is easy to verify that 7 G [a — tt, /3o) (see 
Figure 5.3). 




Fig. 5.3 




Fig. 5.4 



Let a denote the orthogonal distance from to i. Since i necessarily passes between and 
1, it follows that a belongs to the interval (0, 6), where b = — sin 7 denotes the orthogonal 
distance from to the line through 1 with direction e*'''. Let us write f = fr U fi, where 
fr terminates (and // originates) at an inflection point / of /. Since fr G Cr{u,£) and 
fi G Ci{i,v), it follows that 



(5.2) 



\\ff = \\fr 



\' + \\flf> 



|r(w, £) 



\" + \m,v)f 



2 2 

a b — a 
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Let H be the function defined in Lemma 4.2 with Ai = yi, A2 = y2 and = | G (0, 1). 

2 1 

Then (5.2) can be expressed as ||/|| > ^H{ij,), and it follows from Lemma 4.2 that 
H{iJ,) > {1/1+1/2)'^, with equality if and only if fi = yi/{yi +2/2); that is, if and only 
if a = := — sin 7 y^^y^ ■ We therefore conclude that ||/||^ > ^(7) if a 7^ a-^, and 

||/||'>G(7)ifa = aJ' 

Let s^{u,v) denote the set of all non-degenerate right-left s-curves / G S{u,v) with 
inflection direction 7 and let s* (w, v) be the subset of s^{u, v) where a = a^. Let £j be the 
inflection line of a curve in s*{u, v), which is well-defined since all curves in s*{u, v) have 
the same infiection line. We summarize the above discussion, along with the discussion 
preceding (5.1), in the following. 

Theorem 5.7. Assume (5.1). For 7 G F and f G s^(u,f), the following hold. 
(i)\\ff>G{^). 

(n)If\\f\\' = G{-f), thenfes;iu,v). 

The significance of the quantity A(^), when /3 < 0, has been seen in that r{u,ip) is 
directly congruent to \{^)E]^_i^^Qy, we now reveal the significance of A(7). 

Proposition 5.8. Assume (5.1). For ^ G [a — 7r,/3o), r{u,£y) is directly congruent to 
A(7)-E[_f^^o] (^'^d l{£^,v) is directly congruent to A(7)i?[o,t2] . 

Proof. Let Ti and T2 be the similarity transformations such that r{u,£^) — Ti o q] 
and l{£-y,v) — T2 o -^[o.tj]- Since the orthogonal distance from E{—ti) to the real axis is 
^{ti) and the orthogonal distance from to £^ is a^, it follows that the dilation parameter 
of Ti equals = ^ y^^y^ (—sin 7) = A (7). By similar reasoning, the dilation parameter 

of T2 equals |=g = ^(- sin7 + yi sin 7/(^1 + ya)) = A(7). □ 



It was mentioned above that cr(/3) corresponds to the signed distance, in the direction 
e*^, from Pi to 1. The following reveals the significance of (7(7). 

Proposition 5.9. Assume 5.1. For '-y G [o: — 7r,/3o), let Pi and P2 denote the terminal and 
initial points ofr{u,£j) and l{£j,v), respectively. Then (7(7) equals the signed distance, in 
the direction e*'^, from Pi to P2. 

Proof. Let h denote the signed distance in question, and put B = P2 — Pi — he^'^ (see 
Figure 5.4). It follows from Proposition 5.8 that / = r{u,£^) U {l{£^,v) — B) is directly 
congruent to A(7)-E'[_i^ t^]- Since the projected distance, in the direction e*°, from E{—ti) 
to E{t2) equals sinti -|- sint2, it follows that the projected distance, in the direction e*''', 
from the initial point to the terminal point of / equals A(7)(sinti -|- sint2)- Noting that 
the projected distance, in the direction e*''', from to 1 equals cos 7, we deduce that 
A(7)(sinti-|-sint2) = cos 7-/2.. Solving for h and then employing the identity sin^ 5 = sint, 
when 5 = A(£^[o,t]) and t G [0, tt], yields the desired conclusion h = 17(7). □ 

Remark 5.10. Two important consequences of Proposition 5.8 and Proposition 5.9 are: 
1. Let 7 G F, with 7 > a — tt. If 17(7) = 0, then fj := r{u, £-y) Ul{£^, v) has bending energy 
G(-y) and is directly congruent to \{'^)E\^_t^^t^-\. It follows from the latter that fj is of first 
form with to = —ti and t — t2- 
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2. If a {a — tt) > 0, then fa-n '■— fiu,^^) U [Pi,P2] U ^(^7,^) has bending energy G{a — n) 
and is of second form with c = a{a — 7r)/X{a — tt) and t = t2- 

In the following result, we see that (7(7) appears as a factor in the derivative G'{'y). 

Lemma 5.11. Assume (5.1). The function G : F — )■ (0, 00) is continuously differentiable, 

d 1 

has a minimum value Gmin.) cii^d satisfies —GH) = -.o ^il) fof M 7 G F. 
Proof. For 7 e F, we have 

G'ii) = ^{yi + yif - ^(yi + yi){y'Ai) + 2/2(7)) 

sm 7 sin 7 

= ("'^ " ^(-^Vsin(a-7)-^Vsin(^-7))) - ^(^-(t), 

and we note that both A and a are continuous on F and A is positive. If /? < 0, then 
F = [tt — TT, /3] and it is clear that G has a minimum value. On the other hand, if /3 > 0, 
then F = [a — 7r,0), but we note that (^(7) ^ 00 as 7 ^ 0~; hence G has a minimum 
value. □ 

In preparation for the main result of this section, we remind the reader that S{u,v) has 
been partitioned as 

(5.3) Siu,v) = Sl{u,v)U [j s^iu,v), 

where S'i^{u,v) is non-empty only when /3 < 0. 

Corollary 5.12. Let (5.1) be in force, and in case P <0, assume a"(/3) > 0. The following 
hold. 

(i) Ifp < 0, then Wff > Gmin for all f e Sl{u,v). 

(a) If^&T and G{^) > Gmin, then \\f\\ > Gmin for all f e s^{u,v). 

(Hi) IfG{a — n) — Gmin, then a {a — it) > and the curve fa-ir, defined in Remark 5.10, is 

the unique curve, modulo equivalence and elongation of u-turns, in Sa-ni'U'jv) with bending 

energy Gmm- 

(iv) If E {a — 7r,/3o) and G^y) = Gmin, then (1(7) = and the curve f^, defined in 
Remark 5.10, is the unique curve (modulo equivalence) in s^{u,v) with bending energy 
C 

(v) If 13 < and G{f3) = Gmin, then cr(/3) = and fp := r{u,£p), which is of first form 
with to ~ —ti and t — 0, is the unique curve (modulo equivalence) in sp{u,v) with bending 
energy Gmin- 

(vi) The minimum bending energy in S{u,v) is Gmin- 

Proof. Items (i) and (ii) are immediate consequences Proposition 5.6 and Theorem 5.7 
(i), respectively. For (iii), assume G{a — tt) = Gmin- If o-{a — tt) < 0, then it follows 
from Theorem 5.11 that G'{a — tt) < 0, which contradicts the assumption that G attains 
its minimimum at a — tt; therefore, a{a — tt) > 0. It now follows from Remark 5.10 
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that /q-tt has bending energy Gmin- Now, suppose / G Sa-Tv{u,v) has bending energy 
G{a — tt) = Gmin- By Theorem 5.7 (ii), / belongs to s*_^(tt,v) and writing f = fiU /2, 
as in the discussion preceding (5.2), it follows from (5.2) that ||/i|| = ||r(tt,£Q,_7r)|| and 
ll/sll^ = ||Z(£Q;-7r, ■i')!!^- It can then be deduced from the results of section 4 that / is 
equivalent to fce-n or can be obtained (equivalently) from /a-vr by elongation of u-turns. 
We have thus proved (iii). Turning next to (iv), assume 7 G (a — n,Po) and G{j) — Gmin- 
Then ^'(7) = and by Theorem 5.11, we have cr'(7) = 0. It now follows from Remark 
5.10 that fy has bending energy Gmin and the previous argument can be applied to show 
that if / G Sj{u,v) has bending energy Gmin, then / is equivalent to (elongation of 
U-turns is ruled out since curves in Sy{u,v) do not have u-turns). This proves (iv). For 
(v), assume /3 < and G{P) = Gmin- If (^(/3) > 0, then G'{P) > 0, by Theorem 5.11, 
which contradicts the assumption that G is minimized at /3. Therefore, a{(3) = and it 
follows that //3 := r{u^lp) belongs to sp{u,v) = Cr{u,v) (ie Pi = 1). From the discussion 
following Definition 5.3, we know that ||r(w, £^)|| = G{P) and that r{u,£j3) is directly 
congruent to X{l3)E[^ti.o]- From these it follows that WfpW = Gmin and that r{u^£p) is of 
first form, with to = and t = 0. Since the turning angle in r{u, ip) has magnitude less 
than TT, it easily follows from the results of section 4 (more precision, please) that 
is the unique curve (modulo equivalence) in sp{u,v) with bending energy Gmin and the 
proof of (v) is complete. We now prove (vi). It follows from (i). Theorem 5.7 and (5.3) 
that 11/11^ > Gmin for all / G S{u,v). Since the function G has a minimum, there exists 
7 G F such that ^(7) = Gmin, and it then follows from items (iii), (iv) and (v) that is 
a curve in S{u, v) with bending energy Gmin- This proves (vi). □ 

The remaining case, < with (t(/3) < 0, will be addressed in the following section. 



6. Optimal s-curves, part II 



The purpose of this section is to prove the following two results. 

Theorem 6.1. Let < ti < t2 < tt satisfy t2 — ti < tt. Then -E'[ti,t2] is the unique curve 
(modulo equivalence) in S{E(ti), E{t2)) with minimal bending energy. 

Theorem 6.2. In the notation of section 5, let a G (0, tt) and /? < satisfy (5.1) and 
suppose (t{(3) < 0. Then there exist —tt < ti < t2 < such that E^t^^t^^ is directly similar 
to a curve f G S{u,v). Moreover, the curve f is the unique curve (modulo equivalence) in 
S{u,v) with minimal bending energy. 

For t G (0, tt], let if; and 9, as shown in Figure 6.1b, be the positive angles made by the 
chord [0, E{t)] and the segment i?[o,t]- With i denoting the tangent line to E at E{t), let 
p{t) denote the orthogonal distance from to £. 
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Lemma 6.3. Fort G (0,7r), e{t) > ^{t). 

Proof. We will first show, by way of contradiction, that 6{t) 7^ '(/'(i). Assume 6{t) = ip{t) for 
some t G (0, tt). Let T{z) = C1Z + C2 be the congruency transformation which interchanges 
E{t) and 0, and set W = T o E^t^Q-^, where -E'[t,o] denotes the reversal of -Ejo,*]. Since 
e{t) = ip{t), it follows that W belongs to Cz(E(0), E{t)). But Theorem 4.4 asserts that E^^^t] 
is the unique curve (modulo equivalence) in Cz(i?(0), E{t)) with minimal bending energy. 
Therefore, since W and -E'[o,t] have the same bending energy, they must be equivalent. 
However, they cannot be equivalent because W begins with a nonzero curvature, namely 
2sint, while -^lo,*] begins with curvature 0. This proves that 6{t) 7^ -(/'(t) for all t G (0, tt). 
While 0{ti) and both equal 7r/2, a simple computation shows that their derivatives 
satisfy —6'{ti) = ^^'{tv) = 1/d, and it follows that 9{t) > "(pit) for t G (0, tt) sufficiently close 
to TT. Since 6 and -0 are continuous, we conclude that 6{t) > ■i/j{t) for alt G (0, tt). □ 

Lemma 6.4. For t G [0, tt), p{t)^{t) < {2d - ^(t))^ 

Proof. The orthogonal distance p{t) can be formulated as the magnitude of the cross 
product E{t) X |^,|||| which yields 



p{t) = det 
We therefore have 



sint ^(t) 
costv 1 + sin t sin t 



= sin^t - C{t) costVl +sin^t, < t < tt. 



p{t)^{t) - {2d - ^{t)f = (sin^ t + 4d)^{t) - 4d^ - (l + cost\/lTsh^) ^{tf 

< (sin^ t + 4d)^{t) - 4d^ =: g{t), 



where the inequality holds since —1 < cost\/l + sin^ t < 1. We note that (/(tt) = and 

(sin'^ t + 4d\ 
3C{t) cost + ^ . It is clear that g'{t) > for t G (0, 7r/2], and for 

v2 — cost / 

t G {tv/2, tt) (where — cost > 0), we have 

g'{t) = (-3cost)sin2t ("-^(t) + —J^I^-i±M==] > {-3cost) sinH ( -C{t) + ^) , 

V —3 COS ry 2 — cost/ \ 0/ 
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as < -costV2 - cost < 1 on (7r/2,7r). Since < ^{t) < d, it follows that g'{t) > 
for all t G (0, tt) and hence g is increasing on [0,7r]. For t G [0, tt), we therefore have 
p{t)C{t) - {2d - < g{t) < g{7r) = 0, which completes the proof. □ 

In the following, we again use the notation S[^{u^v) (resp. S'^i{u,v)) for the set of all 
non-degenerate left-right (resp. right-left) s-curves connecting u to v. 

Lemma 6.5. For t G (0,7r), the following hold: 
(t) Iffe Sl^{E{0),E{t)), then ||/f > \\E[o,t]f ■ 

(ii) If f e S'^i{E{0),E{t)) ends with a left u-turn, then > ||£;[o,t]f • 

Proof. We will employ the notation and results of the previous section, so in order to 
minimize confusion, we will actually prove the following equivalent formulations: 

{i') life Sl^{E{-t),E{0)), then ||/f > \\E[-t,o]f. 

{ii') If / G S'^i{E{-t),E{0)) begins with a right u-turn, then ||/f > \\E[_t^o]f. 

Let T{z) = C\z -f C2 be the similarity transformation determined by T{E{—t)) = and 
T(0) = 1, and note that T brings the configuration {E{—t),E{0)) to the canonical form 
{u,v) (see Figure 6.2), where u = (0, e*"), v = (l,e*^) with a = 0{t), (5 = —tp{t). Since 
< ilj{t) < 9{t) < TV, it follows that (5.1) holds. Noting that r{ujp) = T o E[_tfi], we 
see that a{(3) = and G{(3) = \\r{ujp)f. For {i'), suppose / G Sl^{E{-t), E{0)). Then 
To / G S'i^{u,v), and it follows from Proposition 5.6 that ||To/||^ > ||r(ti, £^) ||^. Since 
r{u,lfi) = To£'[_^ o] 5 we immediately obtain (i'). Now suppose / G 5'^;(£'(—t), £^(0)) begins 
with a right u-turn. Then To/ belongs to the set Sa-7i-(w, v) defined just above Theorem 

2 2 

5.7, and it follows from this theorem that ||T o /|| > G{a — tt). Since G{I3) = ||r(u,£/3)|| 
and r('u, Ip) = To£'[_^ qj, in order to establish {ii'), it suffices to show that G{a—Ti) > G{(5). 
From Definition 5.3, we have G{I3) = ^{t)^ / smi;{t) and G{a - tt) = {2d - ^{t))^ / sm9{t). 
Referring to Figure 6.1b, we see that sm'i/;{t) = $,{t)/\E{t)\ and sin6'(t) = p{t)/\E{t)\. 
Hence 




Proposition 6.6. Let t G (0, tt). Then -E'[o,t] is the unique curve (modulo equivalence) in 
S{E{0),E{t)) having minimal bending energy. 
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Proof. With Corollary 3.4, Theorem 4.4 and Lemma 6.5 in view, it suffices to show that 
11/11^ > ||£'[o,t]||^ whenever / e Sli(E(0), E{t)) does not end with a left u-turn. Let 

f = fr\-\ f I he as stated, where fr is a right c-curvc and fi is a left c-curve (see Figure 6.3). 
Since fr originates on ^(0), it follows that there exists to G ^) such that fi U E^-^^f^^ is 

a left U-turn. Thus / U E[t^to] belongs to S'^i{E{0), E{to)) and ends with a left u-turn. By- 
Lemma 6.5 (ii), 

11/11^ + \\E[t,to] \f = 11/ LJ E[t^to] \f > \\E[o,to] f = \\E[o,t] f + \\E[t,to] 

2 2 

whence we obtain 11/11 > ||i?[o,t]|| . □ 

Remark 6. 7. By symmetry, it follows from Proposition 6.6 that -^[i is the unique curve 
(modulo equivalence) in S{E{t), E(7r)) having minimal bending energy. 

In the context of the previous section. Proposition 6.6 asserts the following. 

Corollary 6.8. Let direction angles a G (0,7r), /3 < satisfy (5.1) and suppose o'{(3) = 0. 
Then Gip/) > G{(5) for all 7 G [a — tt,P); that is, G{j) is uniquely minimized at j — (3. 

Proof of Theorem 6.1. The extreme cases ti = and t2 = tt have been settled in Propo- 
sition 6.6 and Remark 6.7, respectively, so assume < ti < t2 < tt. By symmetry, and 

2 2 

with Corollary 3.4 and Theorem 4.4 in view, it suffices to show that ||/|| > HE'jt^^tj] II 
whenever / belongs to S'^i{E{ti), E{t2)). Let / be as stated, and let 7 G (— 7r,7r) be the 
direction angle of / at an inflection point /. 
Case i.- 7 G [0,7r). 

Then fUE[t2^n] belongs to 5^;(E(ti), ^(tt)), and it follows from Remark 6.7 that \\fUE[t^^^]f > 

ll^[ti,7r]f , which implies ||/f > \\E[t^^t2]f. 
Case 2: 7 G (-tt, 0) 

Since / begins at E{ti) with a direction angle in (0, tt), there exists a point B on /, between 
E{ti) and /, where / has direction angle (see Figure 6.4). Let us write / = /i U /2, 
where /i terminates (and /2 originates) at B. Let £ be the (horizontal) tangent line to / 
at B, and put q — l{£,E{t2)). Since q and -Ejo^tj] ^-re similar and terminate at the same 
unit tangent vector, it follows that q originates at the point of intersection A between £ 
and the line segment [0, E{t2)]. Moreover, since (7 is at a smaller scale than iJjo^tj], we have 

||g||^ > ||i?[o,t2]ll^- Now, it follows from Proposition 6.6 that ||[^, -B] Ll/2||^ > ||g||^, and 
therefore 

11/11' > \\[A,B]Uf2f > \\qf > \\E[0Mf > P[ti,t.]f • 

□ 

Remark 6.9. By symmetry. Theorem 6.1 remains valid when — tt <ti <t2 < 0. 
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angles a> 5 > Q, with a-\- 5 < tx, there exists ti G (— tt, 0) and t2 > ti such that the chord 
[J(ti), J{t2)] intersects J with interior angles a and d at J{ti) and J{t2), respectively. 



Proof. We refer to Figure 6.5. For t G [— 7r,0], let r(t) denote the direction angle of E{t). 
As t ranges from — tt to 0, T(t) decreases continuously from tt to 0, and it follows that there 
exists h G (— 7r,0) such that r(6) = a. For t G [— tt, 6), let Rt denote the ray emanating 
from J{t) with direction angle r(t) — a and note that since the direction angle is positive, 
Rt intersects J at a unique point J{fx{t)), where fx{t) > t. Let oj{t) denote the interior 
angle, at J{fx{t)), made when the chord [J(t), J(/u(t))] intersects J (the interior angle at 
J{t) equals a by construction). It is clear that oj{t) depends continuously on t G [— 7r,6) 
and tends to as t — )■ 6~. We claim that oj{—tv) > d. If tt) > (ie J{fx{t)) lies on 
[0, oo)), then u{—n) = n — a and the claim follows immediately from the assumption that 
a + d < TV. On the other hand, if /^(— tt) < 0, then a = ijji—fxit)) and a;(— tt) = 9{—n{t)); 
hence, by Lemma 6.3, a;(— tt) > a and now the claim follows from the assumption a > 6. 
By the intermediate value property of continuous functions, there exists ti G [— tt, b) such 
that oj{ti) = 5, and the lemma is proved with t2 = fJ-iti). □ 

Proof of Theorem 6.2. Put 5 = —f3 > and note that the hypothesis of Lemma 6.10 follows 
from (5.1), and we obtain the conclusion of the lemma. We claim that t2 < 0. To see this, 
assume to the contrary that t2 > 0. Let T{z) = ciz + C2 be the similarity transformation 
determined by T(J(ti)) = and T(J(t2)) = 1- It follows from the discussion following 
Definition 5.3 that ToiJj^^ q] = and that cr^P) — \ci\ (^2 — 0) > 0, which contradicts 

the assumption that cr(/3) < 0. Therefore, t2 < and we conclude, from Theorem 6.1 and 
Remark 6.9, that f = T o E^t^^t^^ is the unique curve (modulo equivalence) in S{u, v) with 
minimal bending energy. □ 

7. Proof of Theorem 1.1 

Given any sequence of points Pi, Pm ^ we denote by A{Pi, Pm) the family of 
admissible curves through Pi, Pm such that they are s-curves between any two consec- 
utive points Pi and P^+i. 

The main theorem of this section is as follows. 
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Theorem 7.1. Given any sequence of distinct points Pi,...,P^ G there is a curve 
c G A{Pi, Pm), which minimizes the bending energy in the family and between any two 
consecutive points Pi and Pi+i it is an s-curve. 

First we have to show that 

Proposition 7.2. Let Pi, Pm G be distinct points, then A{Pi, Pm) 7^ 0- 

Proof of Prop. 7.2. We have to construct a curve through the points Pi, Pm such that 
it is an an s-curve between any two consecutive points Pi and Pi+i. Actually, we wiU prove 
a slightly stronger statement, namely we will construct a closed curve passing through the 
points Pi, P2, Pm and back to Pi again. The constructed curve will be a curve which 
is between the points Pi and Pi+i, therefore it will satisfy our definition for curves in 
Section 2. 




Let Vi be the unit vector at Pi that is orthogonal to the line bisecting the angle 
Z{Pi_iPiPi^i) and satisfies the inequality = Z{vi, PiPi^i) < 7r/2 (see Figure 7.1). 
The vector Vi is uniquely determined except in the case when the point P^+i is on the 
half-line starting at Pi and passing through P^-i. In this case we have two choices for Vi 
and we choose one of them. For this to make sense for vi and Vm as well we set Pm+i = Pi 
and Po = Pm- 

From the construction it follows that \Z{vi-i, Pi-iPi) \ < 7r/2 and \Z{vi, Pi-iPi) \ < n/2 
for alH = 1, 2, n (again we set vq = Vm and from this it is not hard to see that Vi-i, Vi 
is s- feasible, that is there is an s-curve Q-i connecting Vi-i to Vi for alH = 1, 2, m). 

The union of these curves will be the desired curve. This completes the proof of the 
proposition. 

Let w, V be s-feasible {S{u, v) 7^ 0) unit tangent vectors. Denote by | |w, f | p the minimum 
of the bending energy in S{u,v). By the previous sections (Theorem 5.1) this minimum is 
assumed by a curve in S{u, v). 

We will need the following proposition which shows that the limit of s-feasible vectors 
is also s-feasible assuming that the bending energy is bounded. 

Proposition 7.3. Let P^ 7^ Pv be different points o/M^. Let UmVn be s-feasible unit 
tangent vectors with base-points P^ and Py respectively, such that limun = u and limf„ = 
V. If Wun.v.nW'^ are bounded then u,v is also s-feasible. 

Proof. Without loss of generality we can assume that P^ = G C and P^ = 1 G C. 
Let a, ttn be the angle of u and respectively, measured from the positive x-axis and 
similarly /3, /3„ be the angle of f ,t'n, respectively. 
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It is easy to see that if Un, f n ^ire s- feasible configurations with — >■ ±7r then | |wn, 'f^n I P 
oo. By reversing the directions of the curves we can show the same way that if /3n — > ±7r, 
then llttnj'i'nlP — >■ OO- 

This imphes not only that an, Pn £ ("TTj for all 't- G N but from the assumption that 
the bending energies of Uttn^'J^nlP were bounded it follows also that a, /3 e (— 7r,7r). 

If q; = then /3 G (— 7r,7r) and it is an s-feasible configuration. 

If a 7^ we may assume, without loss of generality, that a > and > for all 
n e N. It is not hard to see that Un, Vn are s-feasible if and only if G [a^ — tt, tt). Taking 
into consideration that (3 G (— 7r,7r), we conclude that (3 E [a — tt, tt), therefore the limit 
configuration is s-feasible. This concludes the proof of the proposition. 

Next, we show that the bending energy ||tt, v|p is continuous in u and v. For this we 
will need some preparation. 

Let / : [0, L] — )■ C be a unit speed right c-curve with initial and terminal directions 
u, V, where = G C and Py = 1 G C. Let u, v be c-feasible unit tangent vectors with 
the same base-points as u and v, respectively. Treating them as complex numbers set 
a = arg{u), a = arg{u), P = arg{v) and ^ = arg{v), where the argument is in the interval 
(— TT, tt). If any one of the angles were to be equal to tt or — tt, the configuration u, v or u,v 
would not be c-feasible. 

Since we assumed that / was a right c-curve we have a > and /3 < 0. We have the 
following. 

Proposition 7.4. With the notations introduced above let us assume that \f{t)\ < M for 
some M > 1 and a, \ f3\ > rj for some rj > 0. Then, for every e > there is a d{e, M, rj) > 
(depending only on e,M,r]) such that if u,v are c-feasible unit tangent vectors with the 
same base-points as u and v, respectively, and |q: — q:|, |/3 — /3| < 5, then there is a c-curve 
ceC{u,v) with ||c||2 < ||/||2 + e. 

Proof of Proposition 1.1^.. We will modify the curve / at both ends to obtain the curve c. 
The modifications are essentially independent from one another. 

Although it is not true that if u, v are c-feasible, then u, v (or u, f ) will also be c-feasible 
(a counterexample can be found easily by considering the case when / is a u-turn) but it 
is easy to see that one of them tt, v or tt, v will be c-feasible. Let us assume that tt, v is 
c-feasible. The other case can be treated similarly. 

First, we will describe how to modify the curve near the base point of u to obtain a new 
curve c, whose angles with its chord will be a. and ^. We can make the bending energy of 
the new curve as close as we want to the bending energy of / if o: is sufficiently close to a. 

We start with the case when < o: < a. We modify our curve as follows (see Figure 
7.2). 
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First, draw a line with angle a through a point on the negative x-axis such that this 
line is disjoint from our curve /, then translate it horizontally towards / until it makes 
first contact. Denote by P the point on the negative x-axis where the line intersects the 
X-axis and by /(to) the point (or one of the points) on the curve where the line makes 
contact with the curve. Let ci be the c-curve that goes from P to /(to) along the tangent 
line then follows the curve / from /(to) to 1. 

It is clear from the construction that ||ci|p < ||/|p. Moreover, it is easy to see that 
there is a 5i(e, r], M) depending only on e, rj and M such that if r]/2 < a — di < a < a, then 

/ = dist{P, 1) < ^7^^- Re-scale the curve ci with center 1 and with scaling factor l/l 
to obtain the curve c. Since the bending energy changes inversely with the scaling factor 
we obtain 

l|c||^ = ||ci||^/< 11/11^^^^^ = ||/|P + ./2. 
In case a > a we modify our curve c in two steps (see Figure 7.3). 




First, we draw the tangent lines 1^ and to / at and 1, respectively and append 
line segments to / at both ends such that the chord of the resulting curve ci is parallel to 
the X-axis. Denote the distance of the chord of the new curve and the x-axis by d. The 
endpoints of the curve ci are denoted hy P & 1^ and by Q G Notice that the bending 
energy of Ci is the same as the bending energy of /. 

Second, replace the line segment [P, 0] by an arc of a circle that is tangent to ci at and 
intersects the segment [P, Q] at an angle a. Let us denote the intersection point on [P, Q] 
by P G [P, Q]. One can construct such an arc easily and we leave the details to the reader. 
Denote the resulting curve by C2- Let us remark that C2 is a c-curve since u, v was c- feasible. 
For the bending energy of the curve C2 we have ||c2|P = ||/|P + | |arc of the circle| p. There 
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are two parameters to be determined: d and the radius of the circle. 

It is easy to see that there is a 1/2 > d{r]) > 0, such that if d = d{r]), then dist{P, Q) < 

f ^ /2 

j' 2_|_g/4 - Once d is fixed one can find a S2{d, r],e) > such that if a < a < a+62 < tt— 77/2, 
len there is an arc of a circle (simply by choosing the radius large enough) with bending 

energy less than e/4. We can also ensure (by choosing 82 sufficiently small) that P lies 

between P and Q. 

At this point we have ||c2|P < ||/|p + e/4 and I = dist{P,Q) < - 

As before, re-scale the curve C2 with scaling factor l/l and translate it to obtain the 
curve c connecting the points and 1. Since the bending energy changes inversely with 
the scaling factor we obtain 



|e||^ = ||c2||^/< (11/11^ + 6/4)^^±i^ = 11/11 Ve/2. 

Applying the same procedure to the other end of the curve concludes the proof of 
Proposition 7.4 

For given a, (3 & (— 7r,7r) let u,v be unit tangent vectors, with base points and 1, 
respectively, such that a = arg{u), (3 = arg{v). It will be useful to introduce the following 
notations: C(q;, /?) := C{u,v), S{a,P) :— S{u,v) and ||q:,/3|P := which we will 

call the bending energy of the pair a, f3. We say a, f3 are s- feasible if and only if u,v are 
s-feasible. 

Proposition 7.5. For any e > there is a //(e) > such that if \a\, |/3| < n{e), then 
\\a,p\\^<e. 

Proof of Proposition 7. 5. All one needs to do is to construct s-curves with small bending 
energy with initial and terminal directions u and v. The construction is easy and we will 
leave the details to the reader. 

In what follows we will rely heavily on the results of the previous sections, mainly 
sections 5 and 6. In those sections it was always assumed that a and {3 were in "canonical" 
arrangement, that is, a > > 0. However, if we modify a and /3 by a small amount 
the new pair 61 and /3 may no longer be in "canonical" arrangement. The following two 
propositions will help dealing with this situation. 

Proposition 7.6. 

fz;||a,/3||2 = 11-/3, 

(ii) \\a,P\\^ = \\-a,-^\\^ 

(iii) \\a,^\\'' = \\^,a\\\ 

Proof of Proposition 7.6. For any curve in S{a,P) if we reflect the curve with respect to 
the X = 1/2 line and reverse its orientation, we obtain a curve in S{—f3, —a) with the same 
bending energy. This means, that there is a bijection between S{a,(3) and S{—(3, —a), 
which preserves the bending energy. This implies (i). Similarly, reflection across the x-axis 
gives a bending energy preserving bijection between S{a, j3) and S'(— ct, — /3), which yields 
(ii). Combining (i) and (ii) we obtain (iii). 
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Set (cK, j3) 



Proposition 7.7. Let us assume that a > |/3| > hut < a < |^|. 

a) if a < P 

{—P, —a) if a < —P 
If \a — a\,\P — P\ < S, then \a — a\,\P — P\ < S. 

Proof of Proposition 7.7. The proof is elementary and we will leave it to the reader. 

Let us indicate how we will use the previous two propositions in the proof of the next 
one. Assume that a > \/3\, and \a — a\,\/3 — P\ < 5 with a > 5 > 0. If o: < |^|, then we 
will replace a, P with a new pair a, $ as in Proposition 7.7. Then the new pair will be in 
"canonical" arrangement {a> \a — a\,\f3 — (3\ < S and from Proposition 7.6 we have 

\\a, PW^ = \\aj\\^ 

Proposition 7.8. With the notations introduced above let us assume that a, (3 are s- 

feasible satisfying a > > 0. Ifw — r]>a>ri, then for every e > there is a S{€,rj) > 
(depending only on e,ri) such that if \a — a\, |/3 — /3| < 5 and a, (3 are s-feasible, then 
WcxJW' < ||a,/3||2 + e. 

Before we start the proof of Proposition 7.8 let us recall some quantities defined in 
Section 5 (Definition 5.3). Assuming that 7 e [a — tt, /?] fl (— cxd, 0) we have 



2/1(0;, 7) = 2 J Vsin rdr, 

2/2(^,7) Vsin rdr, 

G{a,/3,'y) = 



— sm7 

G'niin("' /^) = min{G'(Q;, ^, 7) : 7 G [a - tt, /3] n (/? - tt, 0)}. 

Recall from Section 3 that d = ^ (tt) = | V sin rdr and define the quantity 70 by 
7o = — sin~^ (sin 77(1 — cos?7)^/(16(i^)) . Since tt — a > ry, from the formulas above one can 
verify immediately that 

(7.1) G{a,p,a-7T)<^ <G{a,p,^) if 70 < 7 < 

sm?7 

This implies that if G{a,/3,'y) = G^j^(a, ^), then 7 < /?* = min{/?, 70}. 

It will be convenient to extend the domain of G{a, 13, 7) to include any 7 G [— tt, 0] 
without changing the minimum G.^y^{a, (3) or the value of 7 at which the minimum is 
assumed. Let us define the set to be = {(a, l3) : rj < a < n — ri, < a, a — tt < 
/3}. For (a, /3) G -fC^, /3* = min{^, 70} and 7 G [-tt, 0) we set 

G'(«,/3,/3*) + 7-^* if /3*<7<0 

^(a, /3, 7) = <( ^(a, /3, 7) if a - tt < 7 < ^* . 

G{a, (3,a — 7r) + a — TT — ^ if — 7r<^<a — tt 
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From the remark following inequality (7.1) and from the construction of G it is clear 
that G.^y^{a^0) = Gj^^j^{a,/3) = min{G(Q;, ^, 7) : 7 G [— 7r,0]}. Moreover G and G 
assumes their minimum at the same points, that is Gj^^j^{a, f5) = G{a,P,^) if and only if 

^min("'/^) = G{a,B,-/). 

The quantity (7(7) will be interesting for us only in the case when ^ = j3. Therefore we 
have 

a{a^(3) — cos (3 -\ — — y/ sm.{a — o: — tt < /3 < 0. 

It is easy to see that (j(a, /?) — )■ 1 as ^ — >■ 0. Therefore we can extend the domain of 
a {a, 13) to the region < /3 by setting 

(7(a,^) = l, if ^>0. 

We can summarize the results of the previous sections as follows. 
If a"(a, (3) < 0, then there is a c-curve in S{a,(3) that minimizes the bending energy. 
If a {a, 13) > 0, then there is an s-curve in S{a,j3) that minimizes the bending energy 
and in that case ||q;,/3|P = G-^y^{a^ (3). 

Proof of Proposition 7. 8. 

Case 1. Let us assume that a{a, /3) < 0. Then /3 < and there is a c-curve / e C(q;, 
which is a segment of rectangular elastica minimizing the bending energy in the family 

C{a,P). 

Let us remark that the condition that / is bounded | < M for some M > 1) will be 
fulfilled for M — 10, since the ratio of length over breadth for any segment of rectangular 
elastica is bounded by 10. 

Since a{a, P) is continuous and tt—t} > a > rj > 0, one can see that there is an rji{rj) > 0, 
such that \f3\ > 771(77). Applying Proposition 7.4 to the curve /, it implies that there is a 
5(?7i,?7, e) > such that if \a - a\, \(3 - (3\ < 6, then \ \a., /3\\^ < \\a, + e/2. 

Case 2. Let us assume that a{a, 13) > 0. In this case | |q;, /3| ^ = ^min('-'^' z^) ~ G{a, /3, 7), 
for some not necessarily unique 7 = 7(0, /5). Prom the remark following inequality (7.1) 
we have 7(0;, /3) e [a — tt, /3*], where /3* = min{;5, 70}. 

Prom the definition of G and G one can see that G{a, /3, 7) is continuous, hence uniformly 
continuous on the region Krj x [— 7r,0]. Therefore, there is a 5i(?7, 70) > such that for 
{a, (3), {a', (3') e K^^ we have 

(7.2) |G(a, A 7) - Gia', p', 7')! < |, if l« - \P - P'l \l - l'\ < ^i- 

Let us assume that |q: — q;|, /3| < 5, where 5 < min{5i, ^} is determined later. We can 
further assume, without loss of generality , that a>\^\. Otherwise, replace a, 13 with a, ^ 
defined in Proposition 7.7. Then we have, according to Proposition 7.7, \a — a\^\^ — ^\ < 5 
and we can prove the statement for 5,^. Since Proposition 7.6 implies ||q;,;5|P = ||q;,;5|P 
the statement is proved for a, (3. 

Part (a). If a {a, /3) > define 7 = 7 (a, (3) to be an angle where G(q;, /3, 7) assumes its 
minimum. Since G and G assumes their minimum at the same points (see the remarks 
following the definition of G) we have 

1 ^1 1' = Gmin(«> ^) = G{a, P, 7) = G{a, P, 7) 
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and 

Taking into consideration (7.2) we have 

\\aj\\^ = = G{aJ,^) < G{aJ,^) < G{a,p,^)+^ = + 1. 

This concludes Part (a). 

Part (b). Assume that a{a,P) < 0. Since d < r]/4 we have {a,P), {a, (3) G Kn/2, which 
is a convex set. Therefore the hne segment [(a, /?), (a, /3)] is also a subset of Since a 

is a continuous function of a, j3 there is a pair (cti, e [(a, ^5), (a, ^5)] with cr(Q;i, j3i) = 0. 

Applying the previous argument for ai, Pi instead of a,^ we obtain 

< ||a,/3||2 + £. 

Since cr(Q;i,/3i) = there is a segment of rectangular elastica with = ||Q;i,/3i|p. 

Moreover we have cti > |, therefore we have, as in Case 1, \/3i\ > 771(77/2). By applying 
Proposition 7.4 we obtain 

llaJlP < \\ai,/3i\\^+^, 

a 5 < 5(e/2, 10, min{77/2, 771(77/2)}) as mentioned in Proposition 7.4. 
Combining this with the previous inequality we obtain 

WaJlf < ||a,/3||2 + e. 
This concludes the proof of Proposition 7.8. 

Next we show that: 
Lemma 7.9. The bending energy \ \u,v\\^ is continuous in u and v. 

Proof of Lemma 7.9. Let e > be arbitrary and u, v be s-feasible. Set a, P as before and 
assume that a> \P\ > 0. 

Case 1. If a < /u(e)/2, where //(e) was given in Proposition 7.5, then choose S < 
/i(e)/2. If \a — a], \j3 — l3\ < 5, then a,\^\ < iJ,{e). Applying Proposition 7.5 implies 

\\a, PW^, \ \a, < e, therefore 

\\\aj\\^-\\a,p\\^\<e. 

Case 2. If tt > o: > //(e)/2, then there is a /ii > such that n — 2//i > a > 2//i. Let 
S = min{/xi, 5(e, //i)}, where S{e, ni) comes from Proposition 7.8. If \a — a\, \P — P\ < S, 
then TT — fii > a,a > fii. li a > \P\ > 0, then applying Proposition 7.8 twice to a, j3 and 
to a, ;S we obtain 

< ||q;,;5||^ + e and | |q;, ^| |^ < | |a, ^| |^ + e, 

hence 
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If q: < 1^1, then, as in the previous proof, replace a, P with a,l3 defined in Proposition 
7.7. Then we have, according to Proposition 7.7, |q: — q;|, |^ — ^| < 5 and we can prove 
the statement for a,p. Since Proposition 7.6 implies ||q:,^|P = ||q!,^|P the statement is 
proved for a, 13. 

This completes the proof of Lemma 7.9. 

Proof of Theorem 7.1. Let G A{Pi, Pm) be a sequence of unit speed curves that 
minimizes the bending energy. Let e M be the points where c"^(t^) = Pi and set 
vf — c^{tf) to be the tangent vectors of at Pj. 

Passing on to a subsequence if necessary we can assume without loss of generality that 

is convergent for every i = 1,2, m as n tends to infinity. Set Vi — lim„^.oo vf. 

Since is a minimizing sequence, ||c"|p is bounded. Proposition 7.3 implies that 
Vi, Vi+i are s-feasible for all i — 1,2, m—l. 

For all i = 1,2, m — 1 let Si G S{vi, f^+i) be a unit speed curve with minimal bending 
energy (see Theorem 5.1 ), and let s be the union of these curves. By construction it is a 
curve, and it is a piecewise s-curve, therefore, s G A{Pi, ...,Pm)- 

We will show that it has minimal bending energy in the family A{Pi, Pm)- 

On the one hand it is clear that 

m 

II s f= E II 11^ . 

On the other hand for i = 1, 2, m—l denote by G S{vf, f J^i) a unit speed curve 
with minimal bending energy in the family S{vl^, f^i) (see Theorem 5.1). We know that 

II c*^ l|2<ll 0^1+"^ 1 l|2 

II II ^11 C [Cj , II , 

therefore, 

m—l m — l 

E II 11'^ E iic"[^r,^m] ii'=ii f. 

i=l i=l 

From Lemma 7.9 we conclude that 

m — l m 

hm y II f= y II ||2=|| s f . 

i=l i=l 

Combining this with the previous inequality we get 

II ^ ll'<ll 1^ 
which proves the claim and the theorem as well. 

Remark. It will be shown in a subsequent paper, that if each of the curves Si are of 
form one (see Definition 5.2), then the resulting optimal curve s is C^. 

Remark. Let us denote by AdosediPi, ■■■jPm) the set of closed curves passing through 
the points Pi,...,P^ such that they are s-curves between any two consecutive points. 
Notice, that the proof of Theorem 7.1 will work equally well for closed admissible curves, 
provided there is at least one closed curve in Aciosed{Pi, Pm)- This is exactly what 
we showed in the proof of Proposition 7.2. Therefore we have the following extension of 
Theorem 7.1: 
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Theorem 7.10. Given any sequence of points Pi,...,Pm £ there is a closed curve 
c e AciosediPi, ■■■iPm), which minimizes the bending energy in the family. 
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